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Abstract 



We investigate the light and heavy meson spectra in the context of the instan- 
taneous approximation to the Bethe-Salpeter equation (Salpeter's equation). 
We use a static kernel consisting of a one-gluon-exchange component and a 



confining contribution. Salpeter's equation is known to be formally equivalent 
to a random-phase-approximation equation; as such, it can develop imaginary 
eigenvalues. Thus, our study can not be complete without first discussing the 
stability of Salpeter's equation. The stability analysis limits the form of the 
kernel and reveals that, contrary to the usual assumption, the confining com- 
ponent can not transform as a Lorentz scalar; it must transform as the timelike 
component of a vector. Moreover, the stability analysis sets an upper limit 
on the size of the one-gluon-exchange component; the value for the critical 
coupling is determined through a solution of the "semirelativistic" Coulomb 
problem. These limits place important constraints on the interaction and 
suggest that a more sophisticated model is needed to describe the light and 
heavy quarkonia. 
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I. INTRODUCTION 



In hadron-structure theory one is interested in describing the hadron as a relativistic 
composite system. To date, most basic properties of hadrons cannot yet be derived from 
QCD — the fundamental theory of the strong interactions. Note, however, that QCD sum 
rules can place some constraints regarding quark-distribution amplitudes in mesons and 
baryons |],||]. With the advent of more powerful computing facilities, lattice gauge theory || 
should provide an increasingly useful means of studying hadronic physics. Yet, at the present 
time it does not provide a convenient framework for a systematic study of a large variety of 
hadronic phenomena. Specifically, with the commission of state-of-the-art facilities, such as 
CEBAFQ other nonperturbative techniques will be required which can be used to incorporate 
phenomena at many different length scales within a single theoretical framework. 

To a large extent our current understanding of hadronic structure is based on the non- 
relativistic constituent quark model Meson properties are well reproduced by a phe- 

nomenological potential consisting of the sum of a short-range one-gluon exchange (OGE) 
component and a long-range confining contribution. A quantitative description of meson 
masses, their static properties and decay rates count among the many successes of the 
model. Generally, one would prefer to have a relativistic and manifestly covariant model. 
For example, a covariant formalism will enable one to relate the wave function (or vertex 
function) in different frames. This becomes essential for calculating hadronic form factors 
at finite momentum transfer. 

The starting point for most relativistic studies of the meson spectrum is the covariant 
Bethe-Salpeter equation ||. The Bethe-Salpeter equation can be regarded as the relativistic 
generalization of the Lippmann-Schwinger equation. However, the Bethe-Salpeter equation, 
being covariant, depends on the zeroth component of the relative four-momentum (i.e., the 
relative energy). Aside from the technical difficulties encountered in handling this extra 
degree of freedom, one must decide in the present case how to generalize the essentially 
nonrelativistic quark- ant iquark potential to four dimensions — a nontrivial task to carry out 
correctly. The difficulty in dealing with the relative energy has led to many different ap- 
proximations to the Bethe-Salpeter equation wherein one works within a three-dimensional 
reduction but attempts to retain fundamental physical principles. There is no obviously cor- 
rect method. Thus, one should study different three-dimensional reductions in the hope of 
isolating model-independent results. Here we work within the instantaneous Bethe-Salpeter 
framework (Salpeter's framework). Although retardation effects and manifest covariance are 
lost, one retains relativistic kinematics, the relativistic character of the potential, and the 
Dirac structure of positive- and negative-energy states. 

The use of the instantaneous approximation commonly employed in the literature en- 
tails other problems besides the loss of retardation and manifest covariance; Salpeter's 
bound-state equation is represented by a nonhermitian Hamiltonian. Indeed, it has been 
recognized that Salpeter's equation is identical in structure to a random-phase- 

approximation (RPA) equation familiar from the study of nuclear collective excitations 1TX 



Thus, it can be rewritten as a hermitian eigenvalue equation — but for the square of the 
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energy. This suggests the possibility of imaginary eigenvalues which would signal the on- 
set of an instability For our problem of interest these imaginary solutions are unphysical 
and their appearance can be precluded by limiting the form of the kernel. This is achieved 
through a stability analysis of Salpeter's equation. We should stress that any study based 
on Salpeter's equation is not complete until the stability analysis is performed. The main 
goal of this paper is to present the stability analysis and to examine the implications for the 
meson spectra. 

We have organized the paper as follows. In Sec. [II], Salpeter's equation is presented and 
the method used to solve it is reviewed. In Section |TTT| , we study the interaction kernel for 
the particular Lorentz structures of interest. The stability analysis for the confining part of 
the kernel is reviewed || and the analysis for the OGE component developed. As a result of 
the stability analysis the form of the instantaneous kernel is constrained. We examine the 
consequences of these constraints on the heavy- and light-meson spectroscopy in Section [IV. 
Finally, our concluding remarks are presented in Section |V|. 



II. FORMALISM 



A. Salpeter's Equation 

In the Salpeter formalism fl^fl , the bound-state spectrum is generated as a solution to the 
instantaneous Bethe-Salpeter equation in the ladder approximation. In this approximation, 
the irreducible Bethe-Salpeter kernel is given by 



V(xi,x 2 ) = V(xi,x 2 ) . 



In Ref. 0, the derivation of Salpeter's equation was illustrated using Green's function meth- 
ods, which yielded the following eigenvalue equation for the Salpeter wave function \ E '- 



(xi,y 2 ) = J d 3 z 1 d 3 z 2 G^ l .^(x 1 ,Z2]z 1 ,y2]E)V^ v/ ( y z 1 ,z 2 )x§ v (z 1 ,z 2 ) 
where the Salpeter wave function is defined by 

xf CT (xi,y 2 ) =< ^o\M^i)My2)\^E > , 



(2) 



(3) 



is the free two-body propagator in the instantaneous approximation (IA), \l/o represents 
the vacuum, and ^ e represents the bound state with energy E. Expanding the fermion fields 
in a single-particle basis and using the properties of the free two-body Green's function then 
gives 



xf CT (xi,y 2 y 



kisi;k2S2 



f4 lSl (xi) Vk 2S2 (y 2 



(ki,k 2 ) + 



Vi 



where the Salpeter amplitudes B and D are defined by 

S slS2 (k x ,k 2 ) = (* |6- 1 (ki)d„(k 2 )|* s ) , 
ZWki,k 2 ) ee (* |<(ki)6t(k2)|*s> , 



(4) 



(5) 
(6) 
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and contain all dynamical information about the bound state. 

Salpeter's equations are more conveniently expressed in an angular momentum basis. 
Projecting out the Salpeter amplitudes, expressing Salpeter's equations for B and D in 
the center of momentum frame, and introducing the partial-wave decomposition of the 
amplitudes in terms of total L and S coupled to the total angular momentum J of the 
bound state 

B SlS2 (k) 

= E {\s l] \s 2 \SM s ){LM L -SM s \JM)Y LML {k)B LS j M {k) , (7) 

SM S LM L JM 1 1 



(-) 1 - Sl - S2 £'- Sl - S2 (k) 
,1 1 

<2 Sl '2 



E (\s^-s 2 \SM s )(LM L -SMs\JM)Y LtML {k)D LSJM {k) , (8) 



SM S LM L JM 

one can write Salpeter's equations in an angular momentum basis as 

[+E-2E k }b LSJ (k)= r^E ( (k-LSJ\V ++ \k>-L>S'J)b L , s ,j{k>) 



+ (k; LSJ\V + -\k';L'S' J)d vs ,j{k')} , (9) 



r°° nk ( 

[-E - 2E k ] d LSJ (k) = j( _ £ | <fc ; LSJ\V~ + \k'; L' S' J)b L , s ,j(k') 

+ (k; LSJ\V—\k'; U S' J)d us >j{k')} , (10) 



with b(k) = kB(k) and d(k) = kD(k). For local interactions, such as the ones considered 
here, the matrix elements of the potential are given by (a sum over greek indices is implicitly 
assumed, and a = 1 — a) 



cs 



(k;LSJ\V ++ \k'; L'S'J) = (k; LSJ\V \k'; L'S'J) 

y.P;a'f}' 



J2(-ir + ^f s , LSJ (k)(S\\[V c (k,k')} a e ;a ,J\S)^%s'Ak'), (11) 



(k;LSJ\V + -\k'; L'S'J) = (k; LSJ\V~ + \k'; L'S'J) 
£(-l)^+^. MJ (A^ (12) 



cs 



where 



Ffs;LsAk) = C a/3 ^)E(«0;/30|AO)(£«SJ||[y A K^) A j o ||L5J) , (13) 

A 



C aft (k) = V4*(-l) a 



Ek + M 
2 Eh 



I E k +M " 
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and E). = \Jk 2 + M 2 , with M the constituent quark mass. (The equal-mass case is con- 
sidered here; however, the extension to unequal masses is straightforward.) The quantum 
numbers L, S range only over the values allowed by J w , and correspond to the usual "non- 
relativistic" quantum numbers, while C, S can take on all values allowed by the coupling to 
J, thus reflecting the role of relativity in the calculation. For E > 0, the amplitudes b, d 
satisfy the RPA normalization condition [|TJ| 



E 

LS 



°° dk 



o (2vr) 3 



d 2 LS j(k) 



1. 



(15) 



Eqs. (|S|) and ([T^) are similar to the equations in Ref. [[TDJ] for the two-fermion case, 
form is used for convenience, as the two can be related by charge conjugation. 



This 



B. The RPA Equation 



Salpeter's equation can be cast in the following compact matrix form: 




E 



B 
D 



(16) 



where the matrix elements of the "Hamiltonian" are given by 



(k; LSJ\H ++ \k'; L'S'J) = (k; LSJ\V ++ \k'; L'S'J) + 2E k (2n) 3 5(k - k')5 LU 5 S s> , (17) 
(k;LSJ\H + -\k';L'S'J) = (k; LSJ\V + ~\k'; L'S'J) . (18) 



One recognizes that Salpeter's eigenvalue equation, as given by Eq. (|T6"D, has the same alge- 
braic structure as an RPA equation |7|-[10|. Having identified the algebraic (RPA) structure 
of Salpeter's equation, the same formalism developed by Thouless in his study of nuclear 
collective excitations [11] will be employed. Salpeter's, and in general any RPA-like, equa- 
tion can be rewritten as a Hermitian eigenvalue equation for the square of the energy ]T5|.p~I 



This implies that while the square of the energy is guaranteed to be real, the energy itself 
might not. The appearance of solutions having E 2 <0 signals, in the context of nuclear col- 
lective excitations, an instability of the ground state against the formation of particle-hole 
pairs — a collective mode with imaginary energy can build up indefinitely. Thouless has 
shown that the stability of the nuclear ground state depends on the Hermitian matrix 



(19) 




being positive-definite — all its eigenvalues must be greater than zero pTT] , p~3|jr^1 . This condi- 
tion is equivalent to requiring that both the sum and difference matrices 

H+ = (H ++ + H+~) , (20) 
H~ = (H ++ - H + -) , (21) 



be positive-definite ||14|| . In this form, the stability condition of Salpeter's equation is reduced 
to finding the eigenvalues of the two Hermitian matrices H + and H~ . Thus, the existence of 
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a single negative eigenvalue, of either H + or H~~ , suffices to signal the instability. It is this 
criterion that was employed in Ref. || to examine the Lorentz structure of the confining 
potential; also, it will be used in this work to set limits on the strong coupling for the 
instantaneous OGE kernel in Salpeter's equation. With these limits in hand, one is then 
able to carry out (see Section [TV]) a study of the meson spectra in the framework of Salpeter's 
equation. 



C. Numerical Solution of Salpeter's Equation 

Salpeter's equation is solved via expansion of the Salpeter amplitudes in a suitable basis, 
thus enabling one to treat the instantaneous confining and Coulomb kernels in configuration 
representation and the relativistic kinetic energy operator in momentum representation, 
where they are respectively local. Here and in Ref. 0, one uses the radial eigenfunctions of 
the nonrelativistic harmonic oscillator, R n L, to expand the two amplitudes B^sj and D^sj 
in terms of unknown coefficients 

W>max 

B LSJ (k) = J2 B nLSJ RnL(k) (22) 

n 

Tl"max 

D LSJ (k) = £ D nLSJ R nL (k) , (23) 

n 

up to n=0,...,n max nodes in the basis, for a finite basis. (Since the interaction is spherically 
symmetric, the magnetic quantum number M only denotes a 2 J + 1-degeneracy and plays 
no dynamical role; hence it can be dropped.) This procedure results in a matrix equation 
for the unknown coefficients B n Lsj and D n isj which can be diagonalized using the method 
developed by Ullah and Rowe [13]. Upon diagonalization, one obtains E 2 and the (previ- 



ously) unknown coefficients from which one then can reconstruct the two amplitudes B^gj 
and D LS j, and, ultimately, the Salpeter wave function x E ■ 



III. STABILITY ANALYSIS OF SALPETER'S EQUATION 

The stability analysis is performed by using potentials V(r) having scalar, timelike- 
vector, and vector Lorentz structures 

{lil 2 , for scalar, 
7 1 °72 , for timelike, (24) 
7172/1 , for vector, 

which are the relevant structures for the meson problem. The analysis is concentrated on the 
pseudoscalar (J 7r =0 _ ) channel; with L=S=0, this is the first channel where the instability is 
likely to develop. For this case, Salpeter's equation, for the reduced amplitudes b(k) = kB(k) 
and d(k) = kD(k), takes the following form 

(+E-2E k )b(k) = j Q ^-p|<A;|V-++|A: / >&(A: / ) + (A; | V+- 1 ^> d(A;') ) , (25) 

roo Ah 1 ( -\ 

(-E-2E k )d(k)= J o ^{(k\V + ~\k')b(k') + (k\V ++ \k')d(k>)j . (26) 
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In spite of the simplicity of the angular momentum content of this channel, the matrix 
elements of the potential are complicated by relativistic corrections. We define the angular- 
momentum components of the potential by 



with j c (x) 
one has 



V c (k, k') = (4tt) 2 / dr j c (kr)V(r)j c (k'r) , (27) 
Jo 

xjc{x) being the Ricatti-Bessel function. For scalar and timelike potentials, 



{k\V ++ \k') 
■E k + M 



2E k 
(k\V + -\k') 
' E k + M 
2Eu 



E k , + M- 

2E k r 



E k , + M" 

2Ey 



l + GQ V (k,k') 1 :2( k ( k ,V 1 (k,k') 



C 2 k + ( 2 k ,}v (k,k')±2( k (yV l (k,k') 



(28) 



(29) 



where the upper (lower) sign in the above expressions should be used for scalar (timelike) 
potentials. For vector potentials, one has 



(k\V ++ \k') 
■E k + M 



2E k 
(k\V + -\k') 
■ E k + M 
2E h 



E k i + M " 

2Ey 



Ey + M" 
2Ey 



Ckd' 



Ck + Cfc' 



3 [Cl + Cl'] Wo(k,k' 



CkQ \v (k,k'), 



where the kinematical variable 

Ck 



k 




(30) 



(31) 



(32) 



E k + M 

has been introduced to quantify the importance of relativity. In particular, for scalar and 
timelike potentials, relativistic corrections arising from the mixing of positive and negative 
energy states (as characterized by (k\V + ~\k'}) appear as O relative to the unmixed 

(Breit) case. This contrasts with the behavior for vector potentials, where both matrix 
elements contain O (1) terms, implying that the contribution of negative-energy states for 
vector potentials will impact results more than the scalar and timelike cases. This can be 
seen by taking the nonrelativistic limit (Cfc, Cfc' - ► 0): 

Vo(/c, k') , for scalar and timelike, 

Vo(k, k') , for vector, 
, for scalar and timelike, 

3Vo(A;, k') , for vector. 

The O (1) term in V + ~ in the vector case stems from the additional spacelike (71 • 72) 
contribution relative to the scalar and timelike cases, contrary to the usual assumption that 
71 -72 induces O (jjp} corrections to the nonrelativistic potential. One should also note that 
the 71 • 72 contribution mixes upper and lower components via a spin-spin term, thus giving 
rise to much stronger splittings than either the scalar or timelike cases. 



{k\V + -\k') 



(33) 
(34) 
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A. Stability Analysis for the Confining Kernel 



In this subsection, the stability analysis carried out in Ref. is reviewed for complete- 
ness. For confining potentials, the Fourier transform [Eq. (p7|) 1 is ill-defined. Hence, in 
examining confinement in momentum space in Ref. 0, the following regularization for the 
spatial part of the potential was employed [|H||T6[|: 

V(r 



are vr = a 



d 2 



-//r 



The Fourier transform of the potential is now well behaved and is given by 



V(k - k') 



drf 



k') 2 + r, 2 



(35) 



(36) 



Evidently, one is interested in studying the stability of Salpeter's equation in the limit of 
7] — > 0. The stability analysis required the explicit evaluation of V + and V~ . These were 
computed with the help of Eqs. fl28|) and (|29| ) 



V + (k, k') 
V-(k,k r ) 



(k\V ++ + V A 
(k\V ++ - V" 



\k') = V (k, k') , 
\k')=V (k,k'){;(k,k') 



(37) 
(38) 



where one introduced relativistic "correction" factors, separately, for scalar and timelike 
confinement 



Uk,k') = 

tv{k,k>) = 



M 2 



kk' K(fc,fc 7 ) 
E k E k ,V {k,k>) 
kk' Vi(k, k') 
E k E ¥ ^ E k E k > V (k, k') 



E k E k i 
M 2 



(39) 
(40) 



For both scalar and timelike confining kernels, H + remained positive-definite, while H 
was positive-definite only for timelike confinement, and not for scalar confinement; a scalar 
confining kernel in Salpeter's equation leads to imaginary-energy solutions — irrespective of 
the constituent quark mass ||. Also, if a mixture of scalar and timelike structure for the 
Lorentz structure of the potential is employed 



r = x 7l 7 2 ° + (l-a;)l 1 l 2 
where x denotes the fraction of timelike structure, then 
V+(k,k') = V (k,k') , 
V~(k, k') = V (k, k') 



M 2 +{2x _ iy kk> V l{ k,k>) 



E k E k i 



E k E k ,V (k } k>) 



(41) 

(42) 
(43) 



Again, H + remains positive definite, while in contrast, H~ is positive definite only for 
x > 1/2. Hence, any mix of scalar and timelike Lorentz structures has stable solutions only 
for x in the interval 0.5 < x < 1. This fact will become important in the study of the meson 
spectra. 
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B. Stability Analysis for the Instantaneous OGE Kernel 



One now considers the (short-range) OGE part of the kernel. The stability analysis is 
performed for a pure Coulomb potential, for both timelike and vector Lorentz structures. 
(Although an OGE kernel of vector Lorentz character will be employed for the spectra 
analysis, the timelike results are presented here for comparison.) The spatial part of the 
(instantaneous) Salpeter kernel is 



Voge(t) 



-a. 



(44) 



with the Fourier transform of Voge{ t ) given by 



W(ik-k'l) 



lim 

?7->0 



—Ana. 



(k - k') 2 + r} 2 



(45) 



where a s is the strong coupling (taken here to be independent of the quark-antiquark sepa- 
ration r, or equivalently, of the momentum transfer Q). 

For the timelike case, one has essentially carried out all necessary manipulations in 
Ref . || . Eqs. (j37D and (|38| ) are completely general for any potential with timelike structure, 
and one only has to calculate the necessary multipoles of the potential in momentum space 
to complete the analysis. One looks to the k = k' limit for which the Coulomb singularity 
structure is manifest. For Vq and Vi, one finds (in the rj — > limit):^ 



V (k, k' = k) 
Vi(k, k' = k) 



2 1 

37T a s - 
2 




(46) 
(47) 



with the leading singularities cancelling in forming the ratio 



Vi(k, k' = k) 
v"o V {k, k< = k) 



lim 



lim 



In I -2- 



+ 0( V 2 /k 2 



(48) 



and one sees that £ v (k, k) — > 1 again. However, in contrast to the confining potential, 
the stability of the RPA matrix in present case is not assured, as the "nonrelativistic" 
equation for H + can become unbounded from below for sufficiently large a s . This implies 
the existence of an upper limit on the effective strong coupling that one can use for a timelike 
Coulomb potential in Salpeter's equation, which affects the determination of wave functions 
and spectra within the model. Thus, it becomes necessary to determine an upper limit on 
a s in order to avoid the instability. Before doing this, the vector case will be considered to 
see if the same problem exists. The matrix elements of the potential are somewhat complex; 
however, V~ and V + are quite simple: 



2 The OGE subscript will be discarded, and future references to V will be taken as indicating the 
Coulomb potential, unless otherwise noted. 
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V + = W (k,k'), (49) 
M 2 



V 



EkE^ 



V (k,k'). (50) 



Since the minus sign in Vo leaves V~ positive, one sees that H~ is actually positive-definite, 
so is bounded from below, and one can concentrate on The eigenvalue equation for H + 
is again the "nonrelativistic" Schrddinger problem, but with a strong coupling four times 
the original a s : 

a = Aa s . 

This is a reflection of the strong upper-to-lower coupling induced by the spacelike component 
of the vector Lorentz structure. Hence, for both timelike and vector cases, one has to 
determine an upper limit for a s at which Salpeter's equation becomes unbounded, noting 
that the upper limit for the vector case will be one-fourth that in the timelike case. 

C. Determination of the Ground State Energy at the "Critical Coupling" 

Evidently, one is interested in the spectrum of the "semirelativistic" Coulomb Hamilto- 
nian 

H = 2Vk 2 + M 2 — — ; a>0. (51) 

r 

As we mentioned in Sec. |I[TJ, the existence of a single negative eigenvalue, of either H + 
or H~ , suffices to signal the instability. If the Hamiltonian ( j51|) is unbounded from below, 
or is bounded but has at least one negative eigenvalue, then the RPA instability develops. 



In 1977, Herbst [17j was able to show that: (a) If a > a c = 4/ir, then the semirelativistic 
Coulomb Hamiltonian is unbounded from below; (b) If a < a c , then all eigenvalues are 
greater than or equal to 0. Specifically, Herbst showed that 



a- 2 



E > 2My 1 - (_j ; for a < a c . (52) 

However, for the RPA equations to be stable one must show that all eigenvalues must be 
positive — a result that does not follow from Eq. fl52|) at the critical coupling. Although our 
primary interest in the spectrum of the semirelativistic Coulomb Hamiltonian stems from 
the stability analysis of Salpeter's equation, the semirelativistic Coulomb problem is still of 



considerable theoretical interest [|i8|J19|| . Thus, in this subsection we present a variational 
analysis that, to our knowledge, represents the best estimate of the ground-state energy 
available in the literature. 

To date, no analytic solution for the ground-state energy of the semirelativistic Coulomb 
problem [Eq. (pip] exists. Thus, we determine an upper bound for the ground-state energy 
(Eq) by using a Rayleigh-Ritz variational method. That is, given a trial wave function |\&), 
an upper bound to the ground-state energy is given by the expectation value of H: 



3 In fact, since V is repulsive, H does not even support bound states. 
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E < (H) = (ty\H\V) ; (tf = 1 . 



(53) 



In practice, one attempts to include the relevant physics by including a set of variational 
parameters into |\&). Then, one minimizes Eq. ( [53"D with respect to those parameters to de- 
termine the variational bound and the "optimal" In configuration space, the variational 
wave function is taken to be 



*(r) 

R(r) 



JV fl -e- 
r 



■yr 



1 + Ci(7r) + c 2 (7r) 2 + 



(54) 
(55) 



where the normalization constant is given by 

(2 7 ) 1+2e 



N, 



E ■ 

n,m=0 



^n—m^m 



2" 



T(l + n + 2e) 



1/2 



(56) 



Note that this form of the wave function is appropriate for the L = channel — where the 
instability should first develop. The expansion for the variational wave function is complete 
and its form is motivated by the analytic solution of the Dirac-Coulomb problem [20]. In 
particular, in the weak-coupling limit (a < 1) one recovers the nonrelativistic result by 
choosing e = 1, 7 = Ma/2, and c\ = c% = • • • = 0, i.e., 



E = 2M(1 
R(r) = 



a' 



M 3 a 3 



-Mar/2 



(57) 
(58) 



In contrast, in the strong-coupling limit (a ~ a c ) one expects that the wave function will 
become localized near the origin, as the energy can benefit from the strong Coulomb at- 
traction. Indeed, the asymptotic behavior of the wave function near the origin is known 
analytically |2T|J2^ 



tan 



a 



(59) 



Note that e — > as a approaches the critical coupling. The variational wave function has 
also an analytic representation in momentum space. That is, 



ry/j \ 1 [2 ™ n T(l + n + e) 



sin 



1 + n + e) tan 1 



,7, 



(60) 
(61) 
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One should note that since the constituent quark mass is the only dimensionful parameter 
in the problem, the dimensionless ratio (H)/M is a function of only the coupling constant 
a. In what follows, all expectation values will be written in units of the constituent mass 
M and expressed in terms of the dimensionless parameter a = 7/M. 

For values of a not too close to the critical coupling (a ^1.25 < a c ~ 1.273) the min- 
imization procedure for (H)/M is straightforward and yields the variational energy and 
optimal parameters (a and e) that are displayed in Table |. Note, the table reflects the ap- 
propriate parameters for the case Ci — c% — ■ ■ ■ — 0. In contrast, the minimization procedure 
is highly nontrivial near the critical coupling (and, thus, for e — * 0) as, both, the kinetic and 
potential energy diverge as 1/e (see Appendix). Hence, in order to compute the variational 
energy for small values of e we perform a Laurent expansion of the expectation value of H, 
i.e., 

(J£t = ^-L + h + h 1 e + --- . (62) 

The first term, h-i, must vanish faster than e ~ (a c — a) 1 / 2 as a approaches the critical 
coupling; otherwise one would contradict Herbst findings. Indeed, we can extract the critical 
coupling (a c = 4/7r) by demanding that the coefficient of the singular term in the series 
vanishes. Note, we have shown in the appendix that h_i vanishes as (a c — a) in the c\ = 
c 2 = • • • = limit (see Eq. [ALp . 

In principle, we could find the variational energy in the a = a c limit by minimizing 
h using the trial wave function (|55|). In practice, however, we can manage only a small 
number of variational parameters. Thus, we proceed by, first, minimizing the expectation 
value of H using only one term in the polynomial expansion in Eq. (f)5|), i.e., we set c x = 
C2 = ■ • • = 0. Note, most of the details related to this minimization procedure are presented 
in the appendix. Next, we compute the variational energy by using a two-parameter (a and 
ci) wave function. This procedure is straightforward but tedious. However, it enables us to 
determine the importance of higher-order "corrections" in the polynomial expansion as well 
as the rate of convergence to the ground-state energy. Our results are summarized below: 

(#> _ f 0.968583, for a = 0.7926 and c x = 0.0000; . . 

Ml 0.968514, for a = 0.9359 and c x = 0.1779. ^ 

Moreover, an initial study with three variational parameters (a, c±, and C2) suggests that 
the two-parameter energy is accurate to, at least, one part per million. To our knowledge 
this represents the most accurate value for the ground-state energy of the semirelativistic 
Coulomb problem presented to date. In this way, our small contribution to Herbst work 
reads: If a < a c = 4/7r, then all eigenvalues (in units of M) are greater than or equal to 
0.968514. Note, our results are consistent with those presented by Raynal and collaborators 
in a comprehensive study of the semirelativistic Coulomb problem ||18|| . Their analysis sets 



lower and upper bounds — differing by less than 1% — for the ground-state energy (0.9650 < 
E Q < 0.9686) at a = a c . 

The remaining question to be answered is how, if at all, does the presence of the linear 
confining potential alter this stability analysis? The answer is that only the value of the 
finite piece is changed, at a = a c . An explicit calculation gives for the expectation value of 
the confining potential (with c\ — ci — ■ ■ ■ — 0) 
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Thus, there is no 0(e~ l ) contribution to the expectation value of (H) coming from the 
confining potential; only a positive contribution, a/2aM 2 , remains at e — > 0. Note that, 
because of the additional dimensionful parameter a, the contribution from the confining 
potential, unlike the Coulomb contribution, depends on the value of the constituent mass. 
Hence, for a Coulomb potential of the form of Eq. (|44T) , stability of Salpeter's equation is 
achieved by demanding that 

J4/7T « 1.273 , for timelike; , , 

as -{ 1/tt w 0.318 , for vector, (b5j 

independent of the constituent quark mass. 



IV. MESON SPECTRA 



The heavy quarkonia (cc and bb) and the light quarkonia (uu and ss) are investigated 
in the instantaneous Bethe-Salpeter framework. The analysis is intended to be qualitative 
in nature, with regard to various effects stemming from the instantaneous approximation; 
hence a simple generalization of the Cornell potential |23 is employed. The virtue of this 
approach is its simplicity. The few parameters of the model, constrained by the stability 
analysis, will emerge from the "best" fit to the meson spectra. One wishes to reproduce 
the interesting physics (e.g. hyperfine structure) through relativistic effects, rather than by 
a fine tuning of many parameters. 

For the heavy quarkonia, the mass spectra are reasonably described, but details such as 
the fine and hyperfine structure are less so. For the light quarkonia, the mass spectra are less 
reasonably described: the pion cannot be accurately modelled within this framework without 
losing the remaining spectra. At the very least, a more sophisticated phenomenology will 
be required to accurately describe static properties of the heavy and light quarkonia. 

For the remainder of the section, the following program is carried out: the various models 
to be considered within the Salpeter framework are defined. Then, phenomenological fits to 
the experimental spectra for heavy mesons are carried out, and a quantitative analysis of 
the various approximations in the Salpeter framework and their effect on the spectra and 
splittings is performed. The latter half of the section is concerned with the light mesons. 



A. Fits to Heavy Mesons 



1. Form of the Interaction 



The spatial part of the potential is based on the Cornell potential; that is, confinement 
is parameterized by a linear potential plus a constant, and the asymptotically- free regime 
is parameterized by an instantaneous OGE potential. Since a good description of the mass 
spectra of the heavy quarkonia can be obtained with a linear confining potential plus the 
nonrelativistic reduction of the OGE piece, this is a natural first choice for the Salpeter 
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equation. (One can, of course, generalize the OGE piece to include running coupling effects 
motivated from perturbative QCD.) 

A mixture of scalar and timelike Lorentz structures for the confining kernel is considered, 
as a number of authors [0,^] have suggested that an admixture of scalar and timelike 
confinement is necessary to reproduce the experimental spectra and splittings. The full 
vector structure is not incorporated in the confining kernel, as this leads to an instability 
similar to the scalar case. The full vector structure in the Feynman gauge is used for the 
OGE kernel. Thus, the instantaneous Bethe-Salpeter kernel is parameterized as 



^(r)r lj2 = [ar + c a ] x ct7i °7 2 + (1 - xjUU 





~a s ~ 




r 



(66) 



where a is the confining string tension, x a controls the mixing between scalar and timelike 
confinement, and a s is the (scale-independent) strong coupling constant. One should note 
that a, c CT , x a , a s , and the quark mass are the only "free" parameters in the approach. From 
the results of the stability analysis, x a and a s are restricted to the following values in the 
Salpeter model in order to have real eigenvalues: 



0.5 < x„ < 1.0 



a s < - « 0.318 

7T 



Recall that a typical value for a s is 0.24 [|2"6" |. One should note that the Breit (V 



= 0) 

model has no such restrictions in principle, saving the possibility of the spectrum becom- 
ing unbound from below if the strong coupling becomes too large, similar to the problem 
encountered for the relativistic Coulomb problem of Section |TTT| . However, for comparative 
purposes, x a is restricted to the same range in both models. 



2. Data and Procedures 

One solves for eigenenergies of the Breit (no coupling between positive and negative 
energy states) and full Salpeter equations, expanding the Salpeter partial-wave amplitudes 
in an oscillator basis of 20 states per partial wave (n max =20, or up to 19 nodes in the 
amplitudes) to insure adequate convergence of the solutions, and with an oscillator parameter 
value (3=0.6 GeV suitable for the heavy mesons. Initially, there are six free parameters in 
the model: a, a s , M c , M b , x a , and c a . The same a, x a , and a s are used to fit both charm 
and beauty. c a is set to zero, assuming the long-range part of the kernel to be less important 
for the heavy mesons. The scalar-timelike mixing parameter x a will be allowed to take only 
the values 0.5 and 1.0, corresponding to equally-mixed scalar and timelike structure and 
pure timelike structure, respectively. Hence, the model for the heavy mesons will only have 
four free parameters, which are determined through minimization of the \ 2 function 

X V, M c , M b ) = £ {E f - E Y pf (67) 

j=l ^ith ~r ®iexp 



using a nonlinear optimization routine [E7J . The Ei exp and <Ji exp are the experimental masses 



and associated errors chosen for the fit, here the first two observed 1 states of b-quarkonium, 
the first two observed 0~ states and the first five observed 1~ states of charmonium. These 
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states are used to determine a set of parameters for each choice of Lorentz structure of 
the confinement. Then, one observes how each of the corresponding spectra agrees with 
the experimental data overall. The E it h and a^h are the calculated masses and their errors 
within the model (which were taken to be 5-10 MeV). The incorporation of a theoretical 
error allows one, in principle, to "force" a better fit to some states, at the risk of possibly 
degrading the fit with respect to the rest of the spectra; however, all states were weighted 
equally in this regard. 

The minimization results for the four parameters for the four models (see Table [H] for the 
appropriate definitions) are summarized in Table [TTT], while the charm and beauty quarkonia 
spectra calculated for each model are summarized in Table [IV] and Table |VJ The experi- 
mental masses were taken from the Review of Particle Properties p8[ (with the exception 
of the h c , which was taken from reference [p^]). Graphical depictions of the spectra for each 
J 71 " channel are shown for the Breit-Timelike model in Figures [l] and |2|. J n states are named 
by their main 2S+1 Lj components. One notes that for all cases, the parameters fall within 
commonly accepted ranges for quark potential models. 

Figure |3] shows the convergence of the ground and first-excited state energies as the 
number of states in the basis is increased; adequate convergence is achieved with 20 states 
in the basis. The numbers quoted for the higher-lying states in Table [TV] and Table [VJ 
should be noted with caution, as typically a 10-20 MeV shift in the energies for the fourth- 
to-fifth eigenstates in going from 18 to 20 basis states is encountered. In particular, the 
larger differences in the table with respect to the experimental data for the higher excited 
states are indicative of an insufficient number of states in the basis. Also, all the models 
appear to be somewhat deficient in comparison to results given by Long J7[ and Spence and 



Vary [26|, which examine both Breit and Salpeter equations for a scalar confining kernel and 
a vector Coulomb kernel. However, those studies allowed more freedom in determining the 
Salpeter solutions. Long utilizes an oscillator basis but minimizes each eigenstate separately; 
each state has a different value of (3 characterizing it, rather than one value for all states. 
While this procedure does minimize the eigenenergies with fewer states in the basis, the 
disadvantage is that the eigenstates are not orthogonal, which would be a problem in the 
calculation of matrix elements, and in ensuring the proper normalization of bound states. 
Spence and Vary use a spline basis 0] as well as an additional interaction (a so-called "Breit" 
interaction) which makes comparison more difficult. However, in their work, solutions with 
imaginary roots (for the light mesons in particular) were discarded on the claim that the 
imaginary roots appear far from the real roots of interest in the complex plane when the 
spline basis is chosen as in the study |26| — a procedure that yields (apparently) stable 
solutions of the Salpeter equation for scalar confinement, in contrast to the results presented 
here. This latter approach differs drastically from the viewpoint adopted in this work, 
which is that the onset of imaginary solutions should indicate that a particular interaction 
is physically inappropriate within the model. (It is rather amusing to note that, when 
examining a kernel that can lead to instabilities, one can "tune" the oscillator basis to get 
apparently stable solutions, but that either shifting the value of (3 or increasing the number 
of states in the basis (or both) reveals the instability.) 
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3. Comparison of Approximations to Salpeter's Equation 



Table [VT] lists the first four eigenenergies for the pseudoscalar and vector channels, using 
the Salpeter- Mixed parameters from Table |T|, to illustrate the various relativistic effects in 
Salpeter's equation. (These are also illustrated in Figure f|.) 

The results show a consistent decrease in the energy of a given level as more relativistic 
effects are included in the calculation. In going from the Schrodinger case (nonrelativis- 
tic kinematics) to the spinless Salpeter case (relativistic kinematics), the energy decreases 
simply because the nonrelativistic kinetic energy increases quadratically for large momenta, 
while the relativistic case only increases linearly; since the potential is the same for both 
cases, the states are more "bound" in the relativistic case. For the Breit case with no lower 
components, the relativistic normalization ( ^g M ) °^ ^ ne ^ ree Dirac spinors suppresses the 
overall potential for large momenta with respect to the nonrelativistic CclS6 ; clS the normal- 
ization varies from 1 in the extreme nonrelativistic limit to | in the extreme relativistic limit. 
However, both attractive and repulsive contributions to the potential are suppressed, and 
the energy is still decreased relative to the nonrelativistic case. The inclusion of Z-graphs in 
the Salpeter case always leads to an added attraction, and consequently to energies reduced 
relative to the Breit case. 

With the introduction of the lower components, the energies are decreased still further, 
for both the Breit and Salpeter cases. In particular, the spacelike part of the vector potential 
(71 ■ 72) makes a large contribution. The decrease in going from Breit to Salpeter is realized 
from the fact that the spacelike part connects the large component of a particle spinor to the 
large component of an antiparticle spinor in V + ~ , hence the contribution from the Z-graphs 
is much larger than that of the direct graphs alone. 



4- Fine Structure Analysis 

One can obtain information on the spin dependence, and thus on relativistic effects, of 
the effective potential for the heavy quarkonia by examining the P-wave fine structure. In 
perturbation theory (which is a good approximation in cc and bb), to O (^2) one can assess 
the relative contributions from a Breit reduction of the potential as 

M ( 25+1 P/ = o,i, 2 ) = Mo + a ss < Si • S 2 > +a LS < L • S > +a T < S 12 > , (68) 

where ass, «ls, and ax arise from the spin-spin, spin-orbit, and tensor components of the 
potential, and 

<S S . 2S(S+l)-3 _[+l, forS=l, 

< Sl • S 2 > -|_ f) for g =0 . (69) 

<LS> =l[J{J+l)-L{L + l)-S(S + l)] = { 7;' (70) 
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< S12 > = < 12 



'Si ■ r) (S 2 ■ r) - -Si • S 2 



> (71) 
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+2, for 3 P x , 

-f, for 3 P 2 , 

+0, for^i, 



S 2 L 2 



3(L- S) 



> 



(72) 
(73) 



(the last two expressions for the tensor component applying for the diagonal elements only), 
with M.q the unperturbed mass. To the extent that perturbation theory is valid for heavy- 
quark spectroscopy, the couplings describe fundamental parameters of nature. Solving the 
four equations [Eq. flSE|)] with the four unknowns yields the data in Tables |V11| and |V111| . 
Perturbatively, the mass of the 1 P± state should be equal to the center of gravity (COG) of 
the 3 Pj multiplet 



1 r 



COG( 3 Pj) = - [5M ( 3 P 2 ) + 3M ( 3 Pi) + M ( 3 P ) 



(74) 



with corrections up to O (i^zj- (One notes that M.q is equal to the COG in the limit of 
a zero-range spin-spin interaction). The spin-spin contribution is, except for the Salpeter- 
Mixed model, an order of magnitude smaller than the tensor and spin-orbit contributions; 
this can be understood by remembering that the spin-spin term in the Breit reduction is a 
contact interaction; since the P-wave states have no support at the origin, spin-spin effects 
are minimized in this channel. (They are not zero here because relativistic corrections in the 
Breit and Salpeter models regularize the contact term.) The X P\ state is off in all models, 
but the error is about 1% at most (0.75% for the Salpeter- Mixed model). For both Breit 
and Salpeter models, an appropriate mixture of scalar and timelike confinement would be 
required for a closer match with experiment for the 1 Pi bb state. 



B. Fits to Light Mesons 

1. Form of the Interaction 

The interaction for the light mesons is taken to be the same form as for the heavy 
mesons. The flavor-independent OGE kernel, however, leads to degenerate 7r and r\ masses; 
one would need to take higher-order diagrams into account that would lead to a flavor- 
dependent interaction, such as annihilation diagrams. These, however, are nontrivial to 
consider in the instantaneous framework and are not treated in this work. It should be noted 
that there are other QCD-based candidates for flavor- dependent qq interactions that have 
been computed by t'Hooft and others from instanton effects |3l|j3^] . Such an interaction has 



been employed in a study similar to this one by Resag et al. |33| for an effective description 
of the light meson spectra. 

One change from the heavy quarkonia is that the constant in the confining kernel is 
permitted as a free parameter. That this is necessary is evinced by Figure |5], which il- 
lustrates the ii-p ground-state splitting, using an up mass M u =M d =0.154 GeV, a string 
tension o~=0.2867 GeV 2 , and initially a strong coupling a s =0.2427 , taken from Spence and 
Vary ||26|| . This coupling (a s ) is then increased to a s ~0.318, the maximum allowed value by 
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the stability analysis in Section |T|. The maximum n-p splitting is less than 400 MeV, still 
about 230 MeV less than the experimental value. The simplest prescription for adjusting the 
masses in order to eliminate the difference is to incorporate a constant c a into the confining 
kernel. This "confinement intercept" has been argued for on other grounds: the necessity 
of regularizing the divergence which appears in treating the linear confining kernel in mo- 
mentum space leads to the appearance of a negative constant in the potential p3 |. It has 
also been argued that the constant can be understood as arising from the gluon condensate 
of the nonperturbative vacuum f3~4* |. For heavy systems, its inclusion is not as important, 
but for light systems (and heavy-light systems) which are affected moreso by the long-range 
potential, its inclusion is necessary for even a fair description of spin-averaged mass spectra 
in Schrodinger and relativised Schrodinger (i.e., spinless Salpeter) approaches; hence, its 
inclusion in the Breit and Salpeter models here is perhaps justified. 

Additional complications arise from the fact that in Salpeter's equation — and in general, 
any two-body quasi-potential equation — a constant term in the confining kernel does not 
solely yield an additive shift in the meson mass spectrum but provides a dynamical con- 
tribution as well, unlike the Schrodinger or spinless Salpeter cases, where one can view it 
as a "negative mass" added to the Hamiltonian. Even with other (repulsive) interactions 
present, c a can be increased to the point where the Salpeter solutions exhibit RPA-type 
instability; the corresponding effect on the Breit solutions is that they become unbounded 
from below, as with the relativistic Coulomb problem of Section |HI| . 



2. Data and Procedures 

As in the case of heavy quarkonia, eigenenergies of the Breit (no coupling between posi- 
tive and negative energy states) and full Salpeter equations were solved for, expanding the 
Salpeter partial- wave amplitudes in an oscillator basis of 20 states per partial wave to insure 
adequate convergence of the solutions, and with an oscillator parameter value /3=0.3 GeV 
suitable for the light mesons. Initially, there are six free parameters in the model: a, a s , 
M u , M s , x a , and c a , where M s is the strange quark mass. x a will be allowed to take only 
the values 0.5 and 1.0, as before, a, c a , a s , and M u are fixed by fitting to the lowest 1~ 
state, the two lowest 1 + states, and the lowest 2 + state for uu, minimizing the chi-squared 
function 

XV, P., a., M u ) = £ {E % h - E T P? (75) 

i=l a ith ' a iexp 

with the errors chosen as for the heavy mesons. M s was then obtained by taking the 
parameters from the fit, and adjusting it to reproduce the <fi mass. 

The minimization results for the four parameters for the four models (Breit-Timelike, 
Breit-Mixed, Salpeter-Timelike, and Salpeter-Mixed) are summarized in Table |D£|, while the 
light and strange quarkonia spectra calculated for each model are summarized in Table [X] and 
Table |XlJ graphical depictions of the spectra for each J n channel are shown for the Salpeter- 
Mixed model in Figures ^| and [7]. J 71 states are named by their main 25+1 Lj components. 

The Salpeter-Mixed model is the best model in this case. The problem with it, however, 
and with the Salpeter-Timelike model, is that in order to fit the p, c a had to be increased 
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to the point where the 7r became unstable. The constituent mass values fall within accepted 
ranges for the Salpeter-Mixed case and the confinement slope is larger than the empirical 
value <r~0.2GeV 2 obtained from spectroscopy. However, it is still in agreement with the 
lattice result cr=0.33to;23GeV 2 [Q; although this last comparison is not very significant 



because of the large error bars. The confinement offset is comparable to that obtained from 
the prescription c a ~—1^fo f34j] . The Breit cases fit the p, but cannot reproduce the tt at all. 
That some difficulty should be encountered in describing the pion in these models should 
not be unexpected. The mass of the 7r is commonly explained in the framework of broken 
chiral symmetry, where it corresponds to an almost massless Goldstone boson; such models 



incorporating chiral symmetry have been investigated by Gross and Milana [35 



In this case, the necessity for the coupling between positive and negative energy states 
for the light mesons is well illustrated; the Z-graphs provide an additional attraction that 
may be necessary in describing the ir as a deeply bound state of a quark and an antiquark 
(although in the present model the attraction is too strong in this channel). This need for 



the V + ~ component in the Salpeter equation agrees with the results of Gara et al. |24| as 
well, albeit for different reasons. The ss states are well-described, with the exception of the 
+ states. It is known, however, that these scalar states can not be represented as simple qq 



states |3(J. Note that in this case the strange-quark mass was adjusted to reproduce the 
the other states are predictions of the model. 



V. CONCLUSIONS 

We have used Salpeter's equation to study the light- and heavy-meson spectra. This 
study was preceded by a stability analysis of Salpeter's equation that proved essential for 
placing limits on the form of the instantaneous kernel. We stress that because of the RPA 
structure of Salpeter's equation a stability analysis must always be performed — regardless 
of the form of the interaction kernel. 

The two main results that emerged from the stability analysis are: 1) the Lorentz char- 
acter of the confining kernel must be timelike or a mixture of scalar and timelike forms, 
contrary to the usual assumption of pure scalar confinement, and 2) an upper limit of 
a s = l/V was set on the strong coupling constant used in the OGE kernel. This value, 
and the corresponding value for the ground-state energy, were obtained from a variational 
solution to the semirelativistic Coulomb problem. To our knowledge this is the best estimate 
presented to date. Having placed limits on the interaction kernel we proceeded to carry out 
a detailed study of the heavy and light quarkonia. 

Static properties of the heavy and light quarkonia within Salpeter's framework have been 
examined using a generalization of the Cornell potential. For the heavy quarkonia the rel- 
ativistic corrections coming into play in the various models were examined. These models 
included Salpeter and Breit approximations having, either, a timelike or a mixture of scalar 
and timelike Lorentz structures for the confining potential. Recall that the Breit approxi- 
mation is obtained by setting V + ~ to zero. Meson masses were adequately described in all 
the models, with the best results obtained using the Breit model with timelike confinement. 
However, a perturbative study of spin-dependent effects (valid for the heavy quarkonia) 
reveals that the fine structure (P-wave splittings) and hyperfine structure ( 3 Si- 3 Di split- 
ting) cannot be simultaneously described in any of the models by simply varying the mixing 
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of scalar and timelike confinement. For the light quarkonia, the mass spectra, except for 
the pion, are best described by the Salpeter model with mixed scalar-timelike confinement. 
However, none of the models were able to describe the pion, or equivalently, the ir-p split- 
ting. For example, in the Breit model the n-p splitting is a "mere" 180 MeV. This difference 
can be pushed up to about 400 MeV in the Salpeter model at the critical coupling; still this 
value is substantially smaller than the experimental splitting of 630 MeV. The additional 
attraction needed to describe the pion would appear to rule out using the Breit models for a 
description of static meson properties. Whether or not this is also sufficient to rule out the 
use of Salpeter 's equation is not clear at this point. Overall, all of the features of meson spec- 
troscopy could not be simultaneously satisfied using the relatively simple kernel employed 
here. At the very least, a more sophisticated phenomenology is required, especially for the 
light mesons and in particular for the pion. It is likely that some form of chirally-invariant 
model will be needed It seems clear, however, that keeping the couplings between 

positive- and negative-energy states is necessary for any realistic description of at least the 
light spectra, and certainly for a combined heavy-light analysis. Moreover, it should also 
be clear that regardless of the form of the kernel, the stability analysis used here must be 
employed in any study that has Salpeter's equation as the underlying dynamical framework. 
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APPENDIX A: EXTRACTION OF THE COEFFICIENTS TO 0(e) IN THE 

VARIATIONAL ENERGY 

For a variational wave function of the forms given in Eq. (55) and Eq. (|6i~|) , one wishes 



to calculate the matrix elements of the kinetic energy operator and potential operator as 
e — > 0. These matrix elements diverge logarithmically for large momenta and small radii in 
their respective integrands; however, one can examine the behavior of the variational energy 
for small e, and determine the critical value of the Coulomb coupling at which the system 
becomes unbounded. In order to compute the variational energy for small values of e we 
expand (T) and (V) as Laurent series |37| 



m e: 



^- = t -^ + t + t 1 e+--- (Al) 

{ ^ = V -^ + v + v 1 e+--- . (A2) 

One wants the leading coefficients t>_i and the zeroth-order ones to, vq as well. By 
minimizing = to + vq we could find the variational energy in the a = a c limit. First we 
calculate the expectation value of H using only one term in the polynomial expansion in 
Eq. (|55|), i.e., we set c\ = c% = ■ ■ ■ = 0. Starting with the potential, one has 
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on 
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(A3) 



where T(z) is the gamma function | ^7[58 |. Note that Eq. ( |A3[ ) is an exact result; there are 
no terms of higher order in e in the series. The kinetic piece requires a little more care; one 
first rewrites the integrand using a standard trigonometric identity: 
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(A4) 



The second integral in Eq. (|A4j) is convergent, so there is no ambiguity in setting e = 
explicitly. The first integral contains a logarithmic divergence for large k; one adds and 
subtracts this divergence to obtain (with a = 7/M and k/M — > k) 
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(T)d is logarithmically divergent for large momenta, for e = 0, so one evaluates the integral 
for e nonzero, then expands about e = up to 0(e). One then has 
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one must be careful to include all 



One should note that in going from Eq. ( |A8| ) to Eq. 
factors in constructing the series expansion. 

The remaining integrals in (T)c are convergent, and one can set e = explicitly. Using 
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,2\2 



8a 

7T JO 



00 „ fc 2 VFTl- WA; 2 + a 2 
dk — 



(k 2 + a 



2\2 



—T c (a) . 

7T 



(A10) 

(All) 
(A12) 
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Then, in the limit that e — > 0, (H)/M can be written as 



M 



a a 



1 8a, 8a ~ . , „ , , 
- + — ln2 + — T c (a) + 0{e) , 
en ix 



(A13) 



where T c (a) is given in closed form for a > by 
f c (a) = < 



ln2-l 



ln2 + lna -l + ii7fc 



a < 1 
a = 1 
a > 1 



(A14) 



By demanding that the coefficient of the singular term h-i in the series vanishes we can 
extract the critical coupling a c = 4 /it. Furthermore, by minimizing the zeroth-order coeffi- 
cient ho we can derive the variational energy (E/M = 0.968583) and the optimal parameter 
(a = 0.7926). 
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FIG. 1. Charmonium mass spectrum for Breit, with timelike confinement and a vector Coulomb 
contribution. The experimental numbers are in the left-hand column for each spin-parity. The 
spectroscopic notation for coupled states is that of the leading component in the calculation. 
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FIG. 2. Beauty quarkonium mass spectrum for Breit, with timelike confinement and a vector 
Coulomb contribution. The experimental numbers are in the left-hand column for each spin-parity. 
The spectroscopic notation for coupled states is that of the leading component in the calculation. 
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FIG. 3. Convergence of the charmonium ground and first excited states for the pseudoscalar 
channel as the number of states in the oscillator basis is increased, for Salpeter with mixed con- 



finement as in Table [II. The basis parameter is /3=0.6 GeV for all cases. 
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FIG. 4. Effects of various approximations to Salpeter's equation, for the parameters in Ta- 



ble III, for the Salpeter mixed-confinement model, for charmonium. and 1 states are in 



each left-hand and right-hand column, respectively, and the spectroscopic notation quoted is the 
dominant component of the calculation. 
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FIG. 5. Ground state splitting of the ir and p mesons, for mixed (scalar +timelike) and pure 
timelike confinement. An oscillator basis with n max =20 and (5=0.3 GeV was used. The model 



parameters were taken from reference [26|, with the resultant splittings indicated by solid lines. 
The dotted lines are with a s =0.318. 
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FIG. 6. Up mass spectrum for Salpeter, with mixed confinement and a vector Coulomb contri- 
bution. The experimental numbers are in the left-hand column for each spin-parity. The spectro- 
scopic notation for coupled states is that of the leading component in the calculation. 
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FIG. 7. Strange quarkonium mass spectrum for Salpeter, with mixed confinement and a vector 
Coulomb contribution. The experimental numbers are in the left-hand column for each spin-parity. 
The spectroscopic notation for coupled states is that of the leading component in the calculation. 
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TABLES 



TABLE I. Minimization of the variational energy with respect to the parameters e and a for 
fixed a, as a approaches a c . 



a 


<H> 
M 


a 


e 


0.100 


1.9975 


0.0495 


0.9861 


n on a 

0.200 


1.9899 


0.0988 


0.9660 


U.4UU 


1.95oz 


0.1956 


0.8977 


0.600 


1.9007 


0.2930 


0.8041 


0.800 


1.8081 


0.3944 


0.6825 


1.000 


1.6583 


0.5071 


0.5222 


1.200 


1.3639 


0.6563 


0.2703 


1.220 


1.3106 


0.6778 


0.2302 


1.240 


1.2435 


0.7033 


0.1817 


1.260 


1.1461 


0.7374 


0.1144 


1.270 


1.0581 


0.7657 


0.05646 


1.273 


0.9933 


0.7854 


0.01533 


1.2731 


0.9874 


0.7871 


0.01170 


1.2732 


0.9786 


0.7897 


0.006226 



TABLE II. Definitions of the various models used in the calculations; "Mixed" and "Timelike" 
refer to the Lorentz structure of the confining kernel. The OGE kernel is always of vector type. 

Mixed Timelike 
Breit [V + ~=0, x a =0.5] [V+-=0, 2^=1] 

Salpeter [V + ~^0, x CT =0.5] [V + ~^0, x a =l] 



TABLE III. Heavy meson parameters for Breit and Salpeter models with mixed 
(scalar+timelike) and timelike confinement, plus an instantaneous vector OGE contribution. An 
oscillator basis with n max =20 and /3=0.6 GeV was employed for the parameter optimization. 



Parameters 


Breit 
Mixed 


Breit 
Timelike 


Salpeter 
Mixed 


Salpeter 
Timelike 


M c [GeV] 
M b [GeV] 
a [GeV 2 ] 

a s 


1.168 
4.573 
0.2991 
0.2678 


1.379 
4.781 
0.1937 
0.4875 


1.251 
4.623 
0.2570 
0.2872 


1.126 
4.561 
0.2743 
0.2165 
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TABLE IV. Charm quarkonia masses in GeV for Breit and Salpeter models with mixed 
(scalar+timelike) and timelike confinement, plus an instantaneous vector OGE contribution, for 



the parameters as in Table |III|. The calculated states are aligned with the observed states according 
to their spin parity, starting from the lowest mass values. An asterisk on an observed value indi- 
cates a state employed in the fits. An oscillator basis with n max =20 and /3=0.6 GeV was employed 
for all states. 







2S+1 Lj 




Breit 


Breit 


Salpeter 


Salpeter 


Meson 


J* 




Mixed 


Timelike 


Mixed 


Timelike 




QT 




2.979* 


2.984 


2.977 


2.958 


2.968 


J/ij> 


l~ 


3 s, 


3.097* 


3.044 


3.111 


3.099 


3.056 


XcO 


n+ 

U 






o.ouu 




3 ^79 


o.ooo 


Xcl 


1+ 


3 P1 


3.511 


3.421 


3.468 


3.412 


3.424 


h c 


1+ 


'Pi 


(3.526) 


3.440 


3.514 


3.468 


3.455 


Xc2 


2+ 




3.556 


3.467 


3.562 


3.499 


3.513 


Vc 


CT 


So 


(3.594)* 


3.645 


3.597 


3.622 


3.647 


1> 


1" 


3 s 1 


3.685* 


3.688 


3.675 


3.693 


3.705 


1> 


1" 




3.770* 


3.726 


3.755 


3.739 


3.732 


1> 


1" 


3 s 1 


4.040* 


4.167 


4.099 


4.134 


4.218 


1> 


1" 




4.159* 


4.192 


4.152 


4.162 


4.234 


4> 


1~ 


3 s 1 




4.566 


4.459 


4.499 


4.659 




1" 


3 i?i 


4.415 


4.583 


4.499 


4.519 


4.669 



33 



TABLE V. Beauty quarkonia masses in GeV for Breit and Salpeter models with mixed 
(scalar+timelike) and timelike confinement, plus a vector OGE (instantaneous Coulomb) contribu- 



tion, for the parameters as in Table |Dj. The calculated states are aligned with the observed states 
according to their spin parity, starting from the lowest mass values. An asterisk on an observed 
value indicates a state employed in the fits. An oscillator basis with n max =20 and /?=0.6 GeV was 
employed for all states. 







2S+l Lj 




T~> "J- 

Breit 


Breit 


Salpeter 


Salpeter 


Meson 


,r 




Mixed 


Timelike 


Mixed 


Timelike 




(T 


'S 




9.440 


9.377 


9.373 


9.432 


T 


1- 


3 SI 


9.460* 


9.468 


9.459 


9.485 


9.488 


Xhn 


0+ 


u 


9.860 


9.821 


9.853 


9.825 


9.807 


Xm 


1+ 


3 P1 


9.892 


9.843 


9.905 


9.850 


9.827 


A. 01 


1+ 


'Pi 

' 1 




9.850 


9.921 


9.841 


9.826 


Xb2 


2+ 


3 P2 


9.913 


9.858 


9.940 


9.865 


9.845 


Vb 


o- 


'So 




9.995 


9.992 


9.949 


9.962 


T 


1" 


3 s, 


10.023* 


10.013 


10.023 


9.999 


9.994 


T 


1" 


3 D, 




10.112 


10.164 


10.100 


10.077 


Xbo 


0+ 


3 Po 


10.232 


10.245 


10.230 


10.214 


10.211 


Xbi 


1+ 


3 Pi 


10.255 


10.263 


10.266 


10.236 


10.229 


Xb2 


2+ 


3 P 2 


10.268 


10.277 


10.293 


10.249 


10.244 


Vb 


o- 


'So 




10.397 


10.350 


10.332 


10.352 


T 


1- 


3 Si 


10.355 


10.412 


10.371 


10.367 


10.375 


T 


\~ 


3 D, 




10.477 


10.462 


10.434 


10.432 


T 


\~ 


3 s l 


10.580 


10.748 


10.655 


10.675 


10.702 



TABLE VI. Various relativistic effects in Salpeter's equation displayed in the pseudoscalar 
and vector channels for charmonium, using the parameters from the Salpeter-mixed model as in 



Table III. S and D label the dominant L-wave component in the calculated energy. Lower on and 
Lower off refer to the lower components of Dirac spinors being present or not, respectively, in the 
calculation. 







Spinless 


Breit 


Salpeter 


Breit 


Salpeter 


r 


Schrodinger 


Salpeter 


Lower-off 


Lower- off 


Lower-on 


Lower-on 


o- 


4.302 (S) 


4.140 (S) 


3.154 (S) 


3.104 (S) 


3.049 (S) 


2.958 (S) 




6.407 (S) 


5.774 (S) 


3.761 (S) 


3.740 (S) 


3.659 (S) 


3.623 (S) 




8.464 (S) 


7.193 (S) 


4.200 (S) 


4.187 (S) 


4.108 (S) 


4.084 (S) 




10.499 (S) 


8.480 (S) 


4.562 (S) 


4.552 (S) 


4.481 (S) 


4.461 (S) 


\- 


4.302 (S) 


4.140 (S) 


3.154 (S) 


3.149 (S) 


3.107 (S) 


3.099 (S) 




6.407 (S) 


5.774 (S) 


3.761 (S) 


3.758 (S) 


3.699 (S) 


3.693 (S) 




6.442 (D) 


5.880 (D) 


3.817 (D) 


3.816 (D) 


3.747 (D) 


3.739 (D) 




8.464 (S) 


7.193 (S) 


4.200 (S) 


4.198 (S) 


4.141 (S) 


4.134 (S) 
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TABLE VII. Spin-dependent parameters for spin-spin, spin-orbit, and tensor contributions to 
the effective potential for charmonium. "COG" refers to the center of gravity of a given multiplet. 
All values listed are in units of [GeV]. 



Parameter 


Experiment 


Breit 
Mixed 


Breit 
Timelike 


Salpeter 
Mixed 


Salpeter 
Timelike 


COG( 3 Pj) 


3.525 


3.440 


3.504 


3.456 


3.464 


'Pi 


3.526 


3.440 


3.514 


3.468 


3.455 


M 


3.526 


3.440 


3.507 


3.459 


3.462 


ass 


-0.001 


-0.0002 


-0.010 


-0.012 


0.009 




0.035 


0.029 


0.063 


0.043 


0.051 


a T 


0.010 


0.005 


0.013 


-0.001 


0.006 


TABLE VIII. Spin-dependent parameters for spin-spin, spin-orbit, and tensor contributions to 


the effective potential for beauty quarkonium. "COG" refers to the center of gravity of a given 


multiplet. The given 1 P\ experimental 


mass is actually the calculated COG. All values listed are 


in units of [GeV]. 












Parameter 


Experiment 


Breit 
Mixed 


Breit 
Timelike 


Salpeter 
Mixed 


Salpeter 
Timelike 


COG( 3 Pj) 


9.900 


9.849 


9.919 


9.852 


9.833 


'Pi 


9.900 


9.850 


9.921 


9.841 


9.826 


Mo 


9.900 


9.849 


9.919 


9.852 


9.833 


ass 


0.0001 


-0.001 


-0.002 


0.015 


0.009 


OiLS 


0.014 


0.010 


0.023 


0.010 


0.011 


ax 


0.003 


0.002 


0.005 


0.002 


0.002 


TABLE IX. 


Light meson parameters for Breit and 


Salpeter models 


with mixed 


( scalar +timelike) and timelike confinement, plus an instantaneous OGE contribution. An oscillator 


basis with n max = 


=20 and /3=0.3 GeV was employed for the data fitting. 




Parameter 


Breit 
Mixed 


Breit 
Timelike 




Salpeter 
Mixed 


Salpeter 
Timelike 


M u [GeV] 


0.2862 


0.3393 




0.3229 


0.4196 


M s [GeV] 
a [GeV] 2 


0.5500 


0.5720 




0.5610 


0.6240 


0.3841 


0.2576 




0.3744 


0.2574 


c CT [GeV] 


-1.448 


-1.089 




-1.427 


-1.157 


a s 


0.2919 


0.3064 




0.2690 


0.2690 
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TABLE X. Light quarkonia masses in GeV for Breit and Salpeter models for the parameters 



as in Table IX, The calculated states are aligned with the observed states according to their spin 
parity. An asterisk on an observed value indicates a state employed in fitting. An "I" indicates 
imaginary eigenvalues. 



Meson 


,r 


2S+l Lj 




Breit 
Mixed 


Breit 
Timelike 


Salpeter 
Mixed 


Salpeter 
Timelike 


7T 




'S 


0.140 


0.627 


0.642 


I 


I 



r 


i- 


3 S1 


0.768* 


0.771 


0.770 


0.768 


0.769 




0+ 


U 


0.983 


1.066 


0.839 


1.014 


0.787 


h 


i+ 




1.232* 


1.203 


1.153 


1.195 


1.148 


Ol 


i+ 




1.260* 


1.169 


1.033 


1.205 


1.081 


7T' 


o- 


% 


1.300 


1.447 


1.368 


1.370 


1.332 


« 2 


2+ 


3 P 2 


1.318* 


1.320 


1.320 


1.319 


1.317 


P 


1" 


3 Di 


1.47 


1.482 


1.316 


1.512 


1.360 


7 


2" 


^2 




1.597 


1.504 


1.622 


1.532 


7T2 


2" 


^2 


1.670 


1.607 


1.566 


1.627 


1.573 


// 

P 


1" 


3 5x 


1.70 


1.556 


1.451 


1.570 


1.460 


IT 


0" 


1 5'o 


1.77 


2.042 


1.940 


2.006 


1.928 



TABLE XI. Strange quarkonia masses in GeV for Breit and Salpeter models for the parameters 



as in Table IX, The calculated states are aligned with the observed states according to their spin 



parity. An asterisk on an observed value indicates a state employed in fitting. 







2S + l Lj 




Breit 


Breit 


Salpeter 


Salpeter 


Meson 


J 77 


M exp t 


Mixed 


Timelike 


Mixed 


Timelike 


7 


Q- 


'So 




0.928 


0.929 


0.640 


0.741 


/o 


0+ 




0.974 


1.401 


1.198 


1.422 


1.164 




1- 


3 s, 


1.019* 


1.019 


1.020 


1.019 


1.020 


fo 


0+ 


3 Po 


1.400 


1.764 


1.570 


1.784 


1.859 


fl 


1+ 


3 P1 


1.426 


1.487 


1.346 


1.501 


1.369 


h 


1+ 


'Pi 




1.515 


1.431 


1.483 


1.412 




2+ 


3 P2 


1.525 


1.577 


1.587 


1.581 


1.536 


fo 


0+ 


3 Po 


1.587 


2.647 


2.398 


2.636 


2.390 


7 


2" 


3 D 2 




1.920 


1.777 


1.927 


1.788 


7 


2" 


'D 2 




1.928 


1.822 


1.923 


1.817 
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